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Abstract. Simple analysis of the leftmost eigenvalue of Ince's equation 
(a boundary value problem with periodicity) resolves an open issue surrounding 
a stochastic Lyapunov exponent. Numerical verification is also provided. 

Let a > b and cr > 0. Consider the stochastic differential equation (SDE) 

rfXi = ° ^ rft + a J "o^ ) ^* ° ^ 

where Wt is scalar Brownian motion with unit variance and o denotes the use of 
Stratonovich calculus. A measure of the stability of such a system is provided by the 
(almost-sure) Lyapunov exponent 

i= lim - \n\Xt\ 

for Xq 7^ 0. Baxendale [U [2] computed a formula for i and, further, proved an 
associated central limit theorem (CLT) 
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No formula has been known for the variance until now. Our calculation of is 
based on the boundary value problem (BVP) 

y"{x) + csin(2a;)?/'(a;) + (A — /iccos(2x)) y{x) = 0, 

y'{0) = y'{7f) = 

due to Ince f3], where c = —{a — h) / a"^ and yU ^ 0. More precisely, if A(/i) is the 
leftmost eigenvalue of the Ince BVP, given yu, then 

In the next three sections, we discuss how Ince's equation arises and the details of 
computing A'(0) and A"(0). The final section is devoted to numerical verification of 
the preceding. 
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1. Equivalence 
Let 7/0"^ denote the rightmost eigenvalue of the differential operator [T] 

Id Ci — b . ,^ . d U , 2 r-2\ 

o^-^ — 5-sin(2a;)— + — (acos x + bsm x) 

2 dx^ la^ dx 

(which is obtained, in part, by projecting the solution of Baxendale's SDE onto 
the unit circle). By the double angle formula for cosine, the BVP we wish to study is 



y'\x) - ^^-^ sm{2x)y'{x) + 



a2 



-^ + 4 + + cos(2a;)) 



y{x) = 



hence c = —{a — h) j a"^ follows, as does 



27 a + 6 
A = T H ^/i; 



hence 

2 2 



A" = .-. = -^A"(0). 

This argument demonstrates the equivalence of Baxendale's setting (in terms of 7'(0) 
and 7"(0)) and our setting (in terms of A'(0) and A"(0)). 

2. First Derivative 
Let y{x,^) denote the solution of the Ince BVP with A = A(/i) that satisfies the 
initial conditions y{0) = 1, y'{0) = 0. Recall that A(/i) is the leftmost such eigenvalue. 
Define ^ 

zix,iJ,) = —y{x,ii). 
Differentiate Ince's equation with respect to /i, yielding 

z" {x, ;u)+csin(2x)/(x, ;u) + (A'(/i) — ccos(2x)) y(x, /i) + (A(/i) — /iccos(2x)) z{x^ fi) = 0. 

Call this (*). Set = 0, yielding 

z"{x, 0) + csin(2x)/(a;, 0) + A'(0) - ccos(2x) = 

because A(0) = and y{x,0) = 1. Multiply both sides of (*) by exp(— | cos(2x)), 
yielding 

[exp(-f cos(2x))/(x,0)]' = (ccos(2x) - A'(0)) exp(-f cos(2x)). 
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Since ^'(0,0) = 0, 

X 

exp(-f cos(2x))2'(x,0) = j (ccos(20) - A'(0)) exp(-f cos(20))rf0. 



Now z'{x, 0) has period vr, hence ^'(vr, 0) = and therefore 

TT 

j (ccos(20) - A'(0))exp(-f cos(2e))rfe = 0. 



It follows that 

TT 

c / cos(2^)exp(-f cos(20))rf^ 



c ' 



A'(0) = = c 



exp(-f cos(2^))(i^ 

where Jq, /i are modified Bessel functions [3]. This reproduces Baxendale's formula 
for i. 

3. ElGENFUNCTION 

In the next section, we will need a formula for z{x, 0). Note that 

y{x, /i) ~ 1 + /i z{x, 0) 

when /i ^ 0, thus a byproduct of our work is an approximation of the leftmost 
eigenfunction of Ince's equation to first order. 
From the following expressions hold: 

oo 

exp(-f cos(2^)) = /o(-f ) + 2 4(-f ) cos(2fc 6), 

k=l 

oo 

exp(-f cos(2^^)) cos(20) = /i(-f ) + ^ (4-i(-f ) + 4+i(-f )) cos(2A; ^) 

k=l 

hence the integrand (ccos(26') — A'(0)) exp(— | cos(26')) is equal to 

(^fc-i(-f) + ^fc+i(-f)) cos(2fc^) - 2c-if^ V4(-f) cos(2fce) 

tit ^0^-2) 



(^'^-i(-f)-27^4(-f)+4+i(-f)J cos(2A:^). 
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Integrating once gives 



exp(— I cos(2a;))/(x, 0) = - ^ sin(2A;x). 

k=l 



Integrating twice gives 

X 

oo . „ 

z{x,0) = ^J2t exp(f cos(2^))sin(2A;^)d^. 

k=i i 

Although we do not explicitly write out the integrals here, they are elementary and 
can be computed symbolically for arbitrary k. 

4. Second Derivative 

Define 

d 

Difi^erentiatc equation (*) with respect to |Li, yielding 

w'\x, jj,) + csin(2a;)w'(x, /x) + \"{fi)y{x, fi) + 2 (A'(yLi) — ccos(2a;)) z{x, fi) 
+ (A(/i) — iiccos{2x)) w{x, ii) = 0. 

Set A* = 0, yielding 

w"{x, 0) + csin(2a;)M;'(a;, 0) + A"(0) + 2 (A'(0) - ccos(2x)) z{x, 0) = 
because A(0) = and y{x, 0) = 1. We again multiply both sides by exp(— | cos(2a;)): 
[exp(-f cos(2x)X(x,0)]'= [2 (ccos(2x) - A'(0)) 0) - A"(0)] exp(-f cos(2x)). 
Since ^'(0, 0) = 0, 

X 

exp(-|cos(2x))«;'(x, 0) = J [2 {ccos{20) - X'{O))z{0,O) - A"(0)] exp(-f cos(2e))d^. 



Now w'{x, 0) has period tt, hence w'{7r, 0) = and therefore 

TT 

J [2 {ccos{20) - A'(0)) z{9, 0) - A"(0)] exp(-f cos{2e))d9 = 0. 



It follows that 

A"(0) = ^, / (ccos(2^) - A'(0)) z{e, 0) exp(-f cos{29))d9 



where z{9, 0) is defined via an infinite series in the preceding section. 
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5. Numerical Verification 
One way to compute A(/i) is to construct the infinite tridiagonal matrix [5] 



M 



( ro V2q-i 

V2qo ri g_2 

qi r2 g_3 

q2 rs q_i 





v 



where r„ = Am? and g„ = (— n + /i/2)c. The leftmost eigenvalue of M is A(/i). 
Another way to compute A(/x) is to solve the continued fraction equation [5] 



A 
2 



Pi 



P2 



P3 



Pi 



14 -12 -A 14 -22 -A 



14 -32 -A 14 -42 -A 14 -52 -A 



where p„ = (— n + /i/2)(n + 1 + jj,/2)c^. 

High precision estimates of the derivatives of A(yu) at zero are found via 



A'(0) 



A"(0) 



A(/i) /i( 

c— 

^o( 

2 



for «i 0. For example, when a = 1, b = —2 and cr = 10, we obtain 

£ = -0.4887503163943852244580286..., 

= 0.0112485762885419873084837... 
as the CLT parameter values. As another example, 

i = 0.3941998582469360577816389..., 

= 0.3841476218435126147382099... 

when (7 = 1 (and a, b remain unchanged). The same results were obtained using 
exact formulation in sections 2 and 4, confirming our work. 
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